Magnetically robust topological edge states and flat bands 
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We study thin strips of three dimensional topological insulators in the presence of a spin-splitting 
Zeeman field. We show that under certain conditions the topological edge states at the sides of a 
strip remain robust against a time-reversal symmetry breaking Zeeman field. For a particle-hole 
symmetric system with Zeeman field perpendicular to the strip we strictly proof that the dispersion 
and the spin-orbital structure of the edge states remains unchanged. When the Zeeman field lies 
parallel to the strip, the Dirac spectrum becomes fiat, but remains intact. Above a critical value 
of the Zeeman field a completely flat surface band develops. We present numerical calculations for 
a realistic tight-binding model of Bi2Se3. These calculations show that even though particle-hole 
symmetry is not strictly fulfilled in this system, these special features are still present. 
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Symmetries play an important role in condensed mat- 
ter physics. Recently, a peculiar state of matter, the 
topological insulator, has been first suggested theoret- 
ically hi, Q and subsequently observed experimentally 
[3|-l8|. In these materials time- reversal symmetry guar- 
antees protection of conducting surface states against 
backscattering, vifhich is a requirement for dissipation- 
less spin transport. Usually, vifhen a time-reversal sym- 
metry breaking perturbation is applied, like a magnetic 
field for example, the surface states are not protected 
anymore. Here, we will show that under certain cir- 
cumstances topological surface states may become robust 
even under a time-reversal symmetry breaking perturba- 
tion. Also, a robust flat band may appear. 

Topological insulators are materials, which have an in- 
sulating energy gap in their bulk, but possess topologi- 
cally protected gapless surface states due to significant 
spin-orbit coupling. Due to Kramers theorem a time- 
reversal invariant perturbation cannot open a gap in the 
Dirac spectrum of the surface states in these systems. 
However, a gap can be opened, when a time-reversal 
breaking perturbation is applied. 

Recently, several materials, like Bi2Se3, Bi2Te3, 
Sb2Te3, or PbBi2Te4 have been experimentally identified 
as three dimensional topological insulators [5 8]. The 
topological protection of the surface states and their spin- 
momentum coupled nature makes these systems partic- 
ularly interesting for spintronics. Several devices have 
been proposed theoretically [9l-[l3j|. In particular, the 
combination with ferromagnetic materials or the appli- 
cation of magnetic fields can lead to interesting effects. 

In the present work we study thin strips of three di- 
mensional topological insulators in the presence of a time- 
reversal symmetry breaking Zeeman field. Zeeman fields 
can be introduced into topological insulators either by 
doping with ferromagnetic dopants jist or by proximity 
to a ferromagnetic material. We show that in this case 
certain additional spin-orbital symmetries are present, 
which are respected by the Zeeman field. These sym- 



metries can make the surface states robust even against 
a Zeeman field. We will also show that at larger Zee- 
man fields a phase transition appears into a state with 
a robust fiat band. Flat bands in solids are particularly 
interesting objects, because the group velocity vanishes, 
allowing for immobile wave packets, localized states, and 
a giant effective mass 18|. 

We start from the generic effective two-orbital tight- 
binding model for a three dimensional topological insu- 
lator suggested by Li et al 15, l3| in the presence of a 
Zeeman field: 



i?(k) =eo(k)l4x4 + ^m,(k)r + J2 



2x2 ' 



(1) 

witheo(k) = C+2L)(l-cos fc:r)+2i:)(l-cos fc^), mo(k) = 
M-2B{1 - cosfc^) - 25(1 - cosky), mi(k) = 2Asmk^, 
and m2(k) — 2 A sin ky. Here, the Dirac F matrices are 
represented by F°d.2 = ^ hx2,Tx (8) ax,Tx (8) ay) in 
the basis of bonding and antibonding pz orbital states. 
The Pauli matrices in orbital space are denoted by 
and the ones in spin space by 0"^. The components of 
the Zeeman field in x, y, and z-direction are denoted 
by Vx.y.z, respectively. The parameters A, B, C, D, and 
M have been derived from bandstructure calculations for 
BizSea in Ref.[l3. 

First, we will study a particle- hole symmetric sys- 
tem with D — Q. In Bi2Se3 the parameter D is 
about a factor of three smaller than B and its mod- 
est infiuence will be considered further below. With- 
out loss of generality we can set C — Q. The Hamil- 
tonian ^ then possesses several spin-orbital symme- 
tries in the absence of a Zeeman field some of which 
remain valid for certain directions of the Zeeman field. 
Specifically, let us consider the following symmetry op- 
erators: 6o = (g) tTa;, 6i = Ta; (g) a^, 62 = Ty (g) Uy, 

©3 = (Tz, and parity P = Tz (g I2x2- In the absence 
of a Zeeman field the Hamiltonian ^ fulfils the follow- 
ing symmetry relations: QoH{kx, fcj,)0Q ^ = H{—kx, ky), 



H[ kx, ky), 



H{kx,ky), and 
PH{kx,ky)P~^ = H{~kx, —ky). In the presence of a 
Zeeman field in z-direction the symmetries Oi, 02, Q3, 
and P are still respected, while 80 is broken. In contrast, 
in the presence of a Zeeman field in x-direction the sym- 
metries 8o, 82, and P are still respected, while 81 and 
83 are broken. All of these symmetry operations respect 
8? = P^ = I4x4, thus their eigenvalues are ±1. 

Let us consider now a strip with finite extension in 
y-direction with y € [0,L], but periodic boundary con- 
ditions in a;-direction. In the absence of a Zeeman field 
the system possesses four topological surface states at the 
two edges y = and y = L, two on each side. These edge 
states can be obtained analytically in small ky expansion 
following the method of Ref. [l^ If L is sufficiently larger 
than the localization length of the edgestates, we may 
consider the edges separately letting L 00. The edge 
state wavefunctions are then linear superpositions of ex- 
ponentially decaying functions of the form e^^ with the 
boundary conditions ipiv = 0) = and ?A(y 00) — 0. 

We will first proof that the dispersion of the edge states 
remains unaffected by a Zeeman field in z-direction. Con- 
sider the Hamiltonian 

Ho{k) - TOo(k)r° + 7712 (k)r2 + v; i2x2 ® tr, (2) 

This is the Hamiltonian ([T]) with C ^ D = except for 
the term. Note, that this Hamiltonian anticommutes 
with the symmetry operation 81, i.e. [Ho,Qi]^ = 0. 
Now assume that I'ikJ is a zero-energy eigenstate of Ho- 
lt is easy to show by direct calculation that such a zero- 
energy edge state exists, if M — 2_B(1 — cos kx) + Vz > 0. 
A second one exists, if M — 2_B(1 — cosfc^,) — Vz > 0. 
Physically, this means that an edge state exist for each 
of the two spin-split bulk bands whenever there is a 
band inversion Due to the anticommutation, |\l/fc^) 
either is an eigenstate of 81 or can be constructed as 
such. Now it turns out that 81 = F^. Therefore l^'fc^) 
will also be an eigenstate of the Hamiltonian ([T]) (with 
C = D = 0) and its energy dispersion will thus be given 
by the prefactor of the term times the 81 eigenvalue, 
i.e. E{kx) = zt2Asmkx. This result is apparently in- 
dependent of Vz (and also of M and B, as long as the 
condition for the existence of the edge state is fulfilled). 
As the Hamiltonian also commutes with 83 we can de- 
duce that the four edge states have to be proportional 
to (1,0,0,1), (1,0,0,-1), (0,1,1,0), or (0,1,-1,0), i.e. 
also the spin-orbital structure of the edge states remains 
unaffected by Vz- It is only the spatial structure that 
changes. Note, that this spin-orbital structure has the 
property that the expectation value of the total spin-z 
component vanishes, because the spin in the two orbital 
subbands is just opposite. This is in contrast to the quan- 
tum spin Hall system HgTe, where the edge states consist 

iSli^. In 






of two orbital subbands having the same spin 
the present case, which is relevant for the edge states 



FIG. 1: (Color online) Numerical dispersions of bulk and edge 
states for \4/M = 0.5, 1.0, and 1.5, respectively. When Vz/M 
reaches 1 two of the four edge states disappear. We have used 
B = A = M = 1, and C = D = 0. 



in a thin Bi2Se3 strip, this special feature explains why 
the edge states are robust against a Zeeman field in z- 
direction. 

The parity operator P relates a state at one edge to 
a state with same energy and opposite kx on the other 
edge. From this we can conclude that the states (1, 0, 0, 1) 
and (1, 0, 0, —1) are localized on different edges and have 
the same spatial structure and the states (0, 1, 1,0) and 
(0, 1, —1, 0), too. For Vz — the symmetry 8082 shows 
that (1, 0, 0, 1) and (0, 1, —1,0) are on the same edge with 
same momentum, but opposite energy. The same is true 
for (1,0,0,-1) and (0,1,1,0). These statements also 
hold for finite Vz even though 80 is broken. 

In order to confirm our analytical results we have di- 
agonalized the tight-binding Hamiltonian ([Ij with C — 
D = and AI = B = A= lona. finite lattice of size 
500 X 200 numerically with periodical boundary condi- 
tions in x-direction. The numerical results, shown in 
Fig. [T] confirm the analytical proof. For Vz/M < 1 we 
find four edge states, two on each side. When Vz/M ap- 
proaches 1, the localization length of two of the edge 
states approaches the width of the system, similar to 
what has been found in Ref. I2I At Vz/M = 1 the bulk 
bands touch the Dirac point (see Fig.[T][b)). When Vz/M 
becomes larger than 1, these two edge states disappear. 
The other two edge states survive, though. The system 
thus makes a quantum phase transition into a quantum 
Hall state. Throughout this transition the dispersion and 
the spin-orbital structure of the edge states remains un- 
changed, while the bulk spectrum changes, of course. 

Let us next consider a Zeeman field in x-direction. In 
this case the dispersion of the edge states changes, be- 
cause I2x2 '^cTx does not anticommute with 81 anymore. 
The symmetries 80 and 82 are still respected. Let \^k^) 
be an edge state of H with energy E. Then, 8082 
is an edge state with energy —E which is localized at the 
same edge. If we have just one edge state on each side, 
I'^k^) = ±8082 l^'fe^) must hold and its energy thus has 
to be strictly zero, i.e. a fiat band appears. As shown in 

Fig. O this happens for Vx > ^M(fc^)2 -f 4^2 sm^{kx), 
where M{kx) = M — 25(1 — coskx). This is a topologi- 
cal flat band as has been discussed recently in graphene 
22 1 and topological superconductors with gap nodes . 
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FIG. 2: (Color online) Numerical dispersions of bulk and edge 
states for a Zeeman field in x-direction with V^jM = 0.6, and 
1.6, respectively. When Vx/M exceeds 1 two of the four edge 
states disappear and a flat band appears. Parameters are the 
same as in figure [T] 
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FIG. 3: (Color online) Numerical dispersions of bulk and edge 
states for Zeeman field in ^-direction with Vz/M = 0.5, 1.0 
and 1.5, respectively. When T4/M exceeds 1, two of the four 
edge states disappear. Here, we have chosen realistic param- 
eters for Bi2Se3, as described in the text. 

Such a flat band appears when the bulk gap closes (23 |. 
as is the case here, too. In the present case the flat band 
is generated by a magnetic exchange field, however, and 
thus could be tuned experimentally. We will show below 
that it can be observed realistically in Bi2Se3. 

In Fig. [2] we show the dispersions of bulk and edge 
states above and below Vx/M — 1. For Vx/M < 1 we 
have four edge states. When Vx/M > 1 two of these 
states disappear and a pair of bulk Dirac points appears. 
These points are connected by a flat band of two edge 
states, one on each side of the strip. 

So far we have discussed the particle-hole symmetric 
case C = D = in Eq. ([T]), in which the two phenom- 
ena of stable dispersion and flat bands become particu- 
larly clear. We will now study parameters, which are 
realistic for Bi2Se3, in order to show that these phe- 
nomena remain and how they are modified. Accord- 
ing to Zhang et al. [13] the parameters are given by 
B = 3.3 eV, M = 0.28 eV, A = 0.49 eV, D = 1.14 eV, 
and C = 0.0068 eV, i.e. particle-hole symmetry is some- 
what broken. The parameter M is the size of the bulk 
gap in the absence of a Zeeman field. The parameter 
t = D/B ~ 0.35 is a measure of the particle- hole asym- 
metry. For t > 1 the system would become topologically 
trivial. 

Following the method of Zhou et al. [l^ we can de- 
termine the dispersion and structure of the edge states 
analytically, when a Zeeman field in z-direction is ap- 




FIG. 4: (Color online) (a) Numerical dispersion of bulk and 
edge states for a Zeeman field in x-direction with Vx/M = 2.2. 
When Vx/M exceeds 1, two of the four edge states disappear 
and an almost fiat band appears, (b) Density of States of 
this system for Vx/M = 4. The flat band appears as a peak 
around E = 0.45 eV. Parameters are the same as in figure [S] 

plied. For the dispersions we find E{kx) = t{M ± 
Vz) ± 2 Ay/1 — P sin kx- The spin-orbital structure of 
the edge states turns out to be of the form (w, 0,0, u), 
(— u, 0, 0, — u), (0, u, u,0), and {0, —v,u,0), where u = 
y/1 — and V = y/l + t^. For t = the particle-hole 
symmetric results above are reproduced. Interestingly, 
the dispersions remain robust except for a constant en- 
ergy shift of ±tVz, which is reduced by the particle- hole 
asymmetry. The Zeeman field does not create a gap in 
the dispersion of the edge states. Also, the spin-orbital 
structure of the edge states neither depends on Vz nor 
on kx- Only the parameter t determines the spin-orbital 
structure. 

In Fig. [3] we present numerical dispersions for a fi- 
nite Bi2Se3 strip in the presence of a Zeeman field in 
z-direction using the parameters of Zhang et al, which 
corroborate our analytical results. Fig. ^a.) shows that 
Vz shifts the energies of the edge states, but does not 
change the shape of the dispersions. When Vz/M exceeds 
1 , two of the four edge states disappear and the system 
makes a quantum phase transition into a quantum Hall 
state. 

When the Zeeman field is applied in x-direction, the 
symmetries Oq and P are still respected, while 82 is 
broken. Let |\E'fc^) be an edge state of H with en- 
ergy E. Then, 0o|^'fc^) is an edge state with energy 
E and momentum —kx, which is localized at the same 
edge. If there exists only one state at each edge we have 
E{kx) = E(—kx), i.e. the dispersion is symmetric in kx- 
The state PQq \'^k^) then is an edge state with same en- 
ergy and momentum kx localized at the opposite edge 
and the edge state dispersion is twofold degenerate. 

In Fig. HJa) we show the numerical dispersion of bulk 
and edge states for a Zeeman field in a;-direction with 
Vx/M = 2.2. When Vx/M exceeds 1, again two of the 
four edge states disappear and an almost flat band ap- 
pears. Due to the broken particle-hole symmetry this 
band has a finite dispersion now. For Vx/M = 2.2 this 
band is clearly separated from the bulk bands. It ends at 
the two Dirac points of the bulk bands, as in the particle- 
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hole symmetric case. In Fig. lUJb) we show the density 
of states of the system for 14 /A/ — 4. Here, the flat 
band appears as a peak around E = 0.45 cV. In order 
to realize this situation experimentally one would need 
Zeeman fields in the range 0.28 eV < T4 < 1 eV, which 
are typical exchange splittings in common ferromagnets. 
The criterion for the appearance of this flat band remains 
the same as in the particle-hole symmetric case above, 

i.e. Vx > \J M{kj.y + AA^ s\r?{kx), because its appear- 
ance is related to the closure of the bulk gap, which is 
unchanged, when Z? > 0. The dispersion of the band 
decreases with increasing value of Vx ■ 

In conclusion, we have shown that spin-orbital symme- 
tries in a topological insulator can create special features 
like robust edge states and flat bands. Specifically, we 
have studied a thin Bi2Se3 strip in a Zeeman field. When 
the Zeeman field is applied perpendicular to the strip we 
have shown that the edge states remain robust even un- 
der the time-reversal symmetry breaking Zeeman field. 
Above a critical value the system can be brought into 
a quantum Hall state via a quantum phase transition. 
When the Zeeman field is applied parallel to the strip 
along its long direction, a topological flat band appears 
above a critical value of the Zeeman field. The combi- 
nation of Bi2Se3 thin films with ferromagnetic materials 
can thus create highly interesting features that can be 
tuned by the polarization of the ferromagnet and lead to 
new spintronics devices. 
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